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ABSTRACT 

The observations of the Kepler space telescope revealed that fundamental- mode RRLyrae stars 
may show various radial overtones. The presence of multiple radial modes may allow us to conduct 
nonlinear asteroseismology: comparison of mode amplitudes and frequency shifts between observations 
and models. Here we report the detection of three radial modes in the star RRLyr, the eponym of the 
class, using the Kepler short cadence data: besides the fundamental mode, both the first and the ninth 
overtones can be derived from the data set. RRLyrae shows period doubling, but switches occasionally 
to a state where a pattern of six pulsation cycles repeats instead of two. We found hydrodynamic 
models that show the same three modes and the period-six state, allowing for comparison with the 
observations. 

Subject headings: stars: individual (RR Lyr) — stars: variables: RR Lyrae — stars: oscillations 



1. INTRODUCTION 

The precise an d continuous observa tions of the Kepler 
space telescope (jBorucki et al.1 [201(1 ) revealed a wealth 
of new features in RRLyrae stars. Among the first dis- 
coveries was the detection of pe riod doubling in three 
stars, including R RLyrae itself (jKolenberg et al.l 120101 : 
ISzabo et al.ll2010f) . Additional frequencies in the Fourier 
transform, significant frequency components beside the 
usual RRLyrae pattern (main period, its harmonics, 
modulation sidelobes and the Blazhko-mo dulation fre- 
quen c y) were also detected in se veral stars (jBenko et al.1 
120101 : IGuggenberger et ail I2012H . Out of these new fea- 
tures, the period doubling phenomenon, the alternation 
of a higher and a lower amplitude pulsation cycle, has 
been successfully modeled. Our hydrodynamic calcula- 
tions revealed that a high-order, 9:2 resonance can oc- 
cur between the fundamental mode and the 9th radial 
overtone str ange mode that lead s to the period-doubling 
bifurcation (Kollath ct al. 2 0111). These results were also 
confirmed bv lSmolec et all (|2011l) . 

However, the nature of the additional frequencies is 
not fully understood yet, except for the ninth over- 
tone which is related to the half-integer peaks and 
causes the period doubling. Some of the remaining 
peaks fall into the ranges where the first or second 
radial overtones are expected, but several stars dis- 
play ev en more peaks, possibly indicating non-radia l 
modes (jBenko et al.l 120101 : IGuggenberger et al.l 120121 ). 
Additional modes were also detect ed in stars observed 
by the CoR oT s pace telescope (iChadid et al.l 120101 : 
iPoretti et aLll2010t: IGuggenberger et al.ll2011l ). The stars 
showing first or second overtone signals are unlike the 
usual double-mode RRd pulsators. The additional com- 
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ponents fall into the mmag range so the amplitude ratios 
are extreme. 

Furthermore, almost all stars that show additional 
modes are Blazhko- variables. Only two stars have been 
found in the Kepler sample that show a small second 
overtone wit h only traces of modulation down t o the 
mmag level (jBenko et al.l 120101 : iNemec et al.1 120111 ). On 
the other hand, about half of the modulated stars in 
the Kepler sample show additional peaks in the Fourier 
spectrum, often together with signs of period doubling. 
Period doubling itself is not a regular process: it is mod- 
ulated and sometimes shows signatures of additional bi- 
furcations. The variations d o not necessarily follow the 
Blazhko-modulation either (jSzabo et al.ll2010l ). There- 
fore it is necessary to investigate the origins of those 
variations. 

2. KEPLER PHOTOMETRY OF RRLYRAE 

RRLyr (KIC 7198959) is one of the brightest stars 
observed by Kepler a mong its asteroseismology targets 
(jGilliland et al. lf20Toh . It was observed in long cadence 
mode in Kepler quarte rs 1 and 2 (or Q l and Q2), with 
29.4-minute sampling. iKolenberg et al.l (|2011l ) provided 
a detailed description of the photometry and the nec- 
essary corrections applied. The brightness of RRLyr 
was initially underestimated in the Kepler Input Cata- 
log (KIC) and some flux was lost from the less-than-ideal 
aperture. However, the corrected average brightness re- 
sulted in an assigned automatic aperture that extended 
beyond the edge of the CCD module and so the pipeline 
rejected the star in quarters 3 and 4. A custom aperture 
had to be developed for the star instead and it is being 
used from quarter 5 onwards. RRLyrae was observed 
through the Kepler Guest Observer PrograrrQ in short 
cadence mode (1 minute sampling) both in Q5 and Q6. 

Maxima of the pulsation cycles had to be determined 
for our analysis. For the long cadence data, we reused 
th e maxima de r ived f rom ninth-order polynomials fitted 
in ISzabo et alj (|2010l ) . For the short cadence data we 
calculated spline solutions instead of polynomials. Long 
cadence data are less well suited to determine the sharp 

6 http://keplergo.arc.nasa.gov/ 
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Fig. 1. — Short cadence Q6 data of RRLyrae. Light grey points are the brighter end of the light curve. Black and dark grey (blue 
and orange in the color version) lines connect the maxima of even and odd pulsation cycles respectively. In the case of the period-six 
visualization (shifted by 0.15 mag), every sixth maxima are connected, such as black/blue lines connect the 6fc, 6fc + 2 and 6fc + 4 points 
(circles), while grey/orange lines connect the 6fc + 1, 6fc + 3 and 6fc + 5 points (dots) where ft = 0, 1, 2. . . . The period-six-like sections are 
indicated with the grey boxes: we identified two occurrences between 55385.3-55406.3 and 55436.9-55448.3. Dashed lines mark the P6 set 
which we analyzed separately (see Figure [3]|. The insert shows the calculated luminosity of a model solution close to the 3:4 resonance. 
Note the similarity of the crossings between the individual blue and orange lines or branches. Model parameters (a„ is the eddy viscosity 
parameter) are: M = 0.64 Mq, L = 45 Lq, T £ ff = 6500 K, a v = 0.0175. The insert covers the same time span as the P6 section. 



maxima of RRLyrae-type variations, but we found the 
results satisfactory for further analysis. 

2.1. Analysis of the photometry 

In principle, return maps would be very useful to 
analyze the variations in period doubling detected in 
RRLyrae stars. Return maps are powerful tools to 
examine the overall dynamics of a system by plot- 
ting the successive values of a selected property against 
each other, e.g. the maximum values of the cycles of a 
(multi)periodic signal. On a return map, single-mode 
oscillation is represented by a single dot (no change), 
period doubling by two dots (small and large cycles), 
and double-mode pulsation is either a discrete number 
of points or a continuous loop, depending on whether 
the period ratio is commensurate or not. Return maps 
can provide an easier way to detect these features than 
conventional Fourier transforms as FTs of RR Lyrae vari- 
ations are dominated by the harmonic peaks of the main 
period and may contain many linear combination peaks. 

In the case of Blazhko-variables however, the modu- 
lation will dominate over almost any other feature. So 
instead of creating a return map, we simply plotted the 
maxima of even and odd pulsation cycles respectively 
(blue and orange lines in Figure [T|) in the light curve 
against time. If the period doubling is constant, the two 
lines should never cross each other. 

However, the pulsation of RRLyr is much more com- 
plicated. The period doubling is omniprese nt, contrary 
to th e first implications from the LC data ()Szab6 et al.l 
2010), but it is clearly not constant. The most curi- 
ous phenomenon occurs — coincidentally — in Q6, the 6th 
quarter of Kepler observations, between BJD 55388 and 



55400 (we refer to this section of the light curve as the 
P6 set, see Figured]). An up-down bounce in each arm 
of period doubling is evident, indicating a possible addi- 
tional bifurcation. Only if we connect every sixth max- 
ima instead of every second can we smooth out most of 
the variations between successive maxima. Here the fun- 
damental mode shows "period-six" behavior, i.e. a pat- 
tern of six cycles instead of two seems to be repeating. 
We found one more possible case and an extension to 
the P6 set where the pulsation temporarily shows traces 
of period-six behavior, indicated by the grey boxes on 
Figure [T] Outside the boxes, both visualizations (ev- 
ery second or sixth maxima) are dominated by period 
doubling, confirming the temporal nature of the period- 
six phenomenon. Reanalysis of the Q1-Q2 data suggests 
period-six pattern at the second half of Ql and the be- 
ginning of Q2 but the gap between the quarters hampers 
the detection. 

The period-six behavior can originate from a resonance 
if two modes are involved in the variation. A resonance 
with the first overtone with a period ratio of 3:4 (0.75) 
would be plausible for example, compared t o the pe- 
riod r atio range of double-mode RR Lyrae stars (jPetersenl 
1973). We note in passing that it would be more cor- 
rect to describe this resonance as 6:8 as the presence of 
the period doubling requires six pulsation cycles to re- 
peat the pattern. Another possible mechanism is a pe- 
riodic window in a chaotic regime. A period bifurcation 
ca scade was indeed foll owed up to period-eight solutions 
in Kollath ~et al.l (|201lD . but we have not found models 
yet where the period-doubling instability alone reached 
chaos. It is much more likely that the period-six behavior 
is connected to a resonance between the period-doubled 
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Fig. 2. — Amplitudes of the highest half-integer peaks at every 
(2n + l)/2/o value (blue points). The Fourier-spectrum of the Q5- 
Q6 data, prewithened with the n/o peaks, is shown in grey. There 
is an excess at 9/2/q, at the pos ition of the ninth overtone, similarly 
to Fig. 8 in lSzabo et"aTI (2010ft ■ 

fundamental mode and another overtone (which configu- 
ration itself can be a periodic window of stability between 
three-mode chaotic solutions). 

2.2. The presence of the first radial overtone in 
RR Lyrae 

The analysis of the Q1-Q2 data of RR Lyr did not 
reveal an y additional peaks bey ond the half-integer fre- 
quencies (jKolenberg et alj|201ll ). The sign of the ninth 
overtone, locked in the 9:2 reso nance, was detect ed in 
those data as an excess at 9/2/n (jSzabo et al.H2010l) . and 
it is present in the Q5-Q6 short cadence data set too 
(Figure [2]). 

However, at least two significant peaks are also de- 
tectable in the expected range of the first overtone in the 
Q5-Q6 data, at / = 2.327 and 2.342 dr 1 , with amplitude 
ratios (Aq/A\, the ratio comparing the amplitude of the 
fo peak of the fundamental mode to the f\ peak(s)) be- 
ing 155 and 170, respectively. The presence of multiple 
peaks indicate amplitude variations and quite possibly 
phase variations as well. If analyzed separately, the Q5 
and Q6 data result in somewhat different values, with 
Ag/Ai — 138 and 102 for a sole significant peak. The 
spectrum of the P6 set (lower panel in Figure [3]) shows 
a stronger signal with Aq/A\ = 67. The amplitudes are 
limited to a few mmags in all subsets, but thanks to the 
photometric precision, the individual peaks have signal- 
to-noise ratios of at least 4-6, the lowest values corre- 
sponding to the spectrum of the Q5 data only, where we 
did not detect period-six-like stages. 

Period ratios (or fo/ fi frequency ratios) vary between 
0.753-0.758 for the different subsets and/or peaks: for 
the P6 set it is 0.756. These ratios are close to, but not 
exactly 0.750, the value corresponding to a 3:4 resonance. 
The period ratio difference might arise from the uncer- 
tainty in the frequency determination, either caused by 
the frequency variations connected to the Blazhko effect 
or, in the case of the P6 set, the shortness of the data 
set. But it is more likely that the system temporarily ap- 
proaches the vicinity of the resonance and the frequency 
difference reflects the distance from it. An additional 
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Fig. 3. — Detection of the first overtone (/i) in the Kepler mea- 
surements of RRLyr. The original spectra are in grey and the 
spectra prewhitened with the main frequency (fo) and its harmon- 
ics (n/o) are in black (blue in the color version). The upper panel 
shows the entire Q5— Q6 data while the lower panel shows only the 
P6 set from Q6 (see FigurefTJ. A forest of peaks is evident at 3/2/o 
in the upper panel, indicating the strong variability of the period 
doubling. 

mode w ith the same ratio (0.754) was detected in UZ 
Vir too (jSodor et al.ll2012ft . 

This is the first time that three radial modes have been 
unambiguously detected in an RR Lyrae star: the funda- 
mental, the first overtone and the third being the ninth 
overtone which is responsible for the period doubling. 
The detection is also strongly supported by the nonlin- 
ear models, described below. It is also clear however 
that the inclusion of period doubling and the first over- 
tone still does not describe the light variations of RR 
Lyr perfectly. The two branches of the period doubling 
cross each other and reverse, with intervals ranging from 
a few pulsation periods to about 30 days which is close to 
an entire Blazhko-pcriod (~ 39 days). These variations 
could be explained by yet another, physically relevant 
frequency close to the half-integer peaks (as it modulates 
the period doubling), or coul d be internal variability, a s 
found in BL Herculis models (jSmolec fc Moskalik | [2012). 

3. MODEL CALCULATIONS 

The modulations and possible bifurcations in period 
doubling raise further questions regarding the dynam- 
ics of the system. To study the the possible pro- 
cesses that might be responsible for this dynamical be- 
havior we conducted a broad nonlinear model survey 
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Fig. 4. — First return map of a three-mode model. The figure 
shows maxima of the luminosity variations plotted against each 
other. Here we show the maxima of the fcth versus the (fc + 2)th 
cycles instead of successive ones to separate the smaller and larger 
cycles of period doubling. Although three radial modes are present, 
the ninth overtone is locked in a resonance with the fundamental 
mode, and it is only detectable through the presence of period 
doubling, hence we observe two loops instead of a single one, both 
representing the double-mode pulsation between incommensurate 
fp and /i frequencies. Model parameters are the same as in Figure 

m 

with the Florid a-Budapest turbulent convective hydro- 
dyna mic code (|Kollath fc Buchlerl 120011: IKollath et all 
2002). The results of this survey will be published in 
a subsequent paper (Kollath et al., in prep.), here we 
only summarize the relevant facts. 

We found that period doubling changes the stability 
of the fundamental mode. The mode in itself is stable 
against the perturbations of the first overtone, but unsta- 
ble against the ninth overtone and the period-doubling 
bifurcation occurs. The period-doubled mode is, how- 
ever, unstable (or at least marginally stable) against 
the perturbations of the first overtone. The result is a 
new pulsation state that looks like a new kind of dou- 
ble mode behavior but is in fact created by the inter- 
action of three different radial modes. In some of the 
models, this two-period state bifurcates further to res- 
onant or chaotic solutions. Detailed properties of these 
models, their occurrence along the instability strip will 
be investigated in forthcoming papers. We note that 
the two chao t ic mo del solutions that were analyzed by 
iPlachy et al.1 (|2012l ) turned out to be three-mode solu- 
tions as well. 

The various dynamical states are best to represent with 
return maps. In case of simple period doubling, we will 
have two dots, corresponding to the lower and higher 
amplitude cycles. If two different modes are present, the 
return map will either show a set of points or a loop, 
depending on the phase relations (e.g. are the two modes 
in resonance or not). 

Take for example the two loops in Figure 0] which is the 
same model as in the insert in Figure [T] The presence 
of a loop at lower values and another at higher values 
indicate that both period doubling and an additional pe- 
riodicity is present in this system, a state that resembles 
the Kepler observations of RR Lyrae. 

3.1. Comparison with the models 

The detection of three different radial modes in an RR 
Lyrae star may allow for nonlinear asteroseismic stud- 
ies: comparison of the observations not only with linear 



pulsation periods but also with mode amplitudes, reso- 
nances and frequency shifts. 

Our first calculations focused on the possible 3:4 res- 
onance. The corresponding 0.750 period ratio, and es- 
pecially the exact observed period ratios in RRLyr be- 
tween the first overtone and the fundamental mode are 
a challenge: they lie at the edge of the accessible range 
for the classical doubl e-mode pulsation i n the Petersen- 
diagram (see Fig. 6 in lSzabo et ail (|2004D ). We have not 
succeeded in finding models with an exact 3:4 resonance 
yet, though we found resonant solutions with higher inte- 
gers corresponding to lower period ratios. These models 
can help to identify possible (near-) resonances between 
the fundamental mode and the first overtone in other 
RR Lyrae stars. 

However, models with a linear period ratio of 0.747 in- 
stead of 0.750 do in fact approach the vicinity of the 3:4 
resonance. These calculations indicate that RRLyr only 
approaches the resonance during the period-six phases 
but does not lock into it. Although the models favor so- 
lutions below the 0.750 frequency ratio instead of above 
(as in the case of RRLyr), the topology, like the dis- 
tances and crossings of the six branches are quite similar, 
especially in the P6 section (see insert in Figure Q]) and 
support our findings. Swapping of either three branches 
simultaneously was not observed in the models, confirm- 
ing that yet another mechanism must be behind some of 
the remaining features, like the modulation of the period 
doubling. 

We also note that the difference between linear and 
nonlinear period ratios in the models approaches 0.003 
and sometimes (especially in resonant states) up to 0.005. 
This nonlinear shift is comparable to the observed period 
ratio range of the Petersen-diagram of the normal double- 
mode RR Lyrae (RRd) stars. 

4. SUMMARY 

The detection of so many interesting features in a 
star of long observational history remarkably illustrates 
the advantages o f spac e-based photometry. A decade 
ago, iSmith et al.l (|2003l ) noted that the residuals after 
fitting the primary components of the spectrum (the 
triplet at the fundamental mode, their harmonics and 
the modulation frequency) exceeded the expected obser- 
vational error. The first observati ons of Kepler revealed 
the period-doublin g phenomenon (|Kolenberg et al.lfeOlOl : 
ISzabo et al.l [20Toh and now the first overtone was also 
detected in the short cadence data. The results and im- 
plications of the analysis can be summarized as follows: 

• Three different radial modes were detected in 
RRLyr, the eponym of its class: the funda- 
mental mode and the first and ninth overtone. 
However, the ninth overtone is locked in a reso- 
nance with the fundamental mode and can be de- 
tected only through t he presence of period doubling 
(IKollath et alJl20ll . causing the observed varia- 
tions to resemble double-mode behavior. 

• The first overtone has a very small amplitude (few 
mmags), in contrast with the fundamental mode 
or the classical double-mode RRd stars, explaining 
the non-detection from the ground. 
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• A state resembling period-six bifurcation was de- 
tected in the light variations of RRLyr. This be- 
havior is most likely caused when the system ap- 
proaches the vicinity of a 3:4 resonance between 
the fundamental mode and the hrst overtone. 

• We have been able to reproduce the three-mode 
state in nonlinear hydrodynamic models. Although 
we have not been able to model the 3:4 resonance 
or the higher than 0.750 period ratio of RRLyr, 
the amplitude ratios can be reproduced and higher- 
order resonances have been found as well. Further- 
more, near-resonance models below the 0.750 value 
have very similar characteristics, raising the hopes 
for applying nonlinear asteroseismology to Blazhko 
RR Lyrae stars. 

• Even the inclusion of three different modes can- 
not explain all the variations we see in the period 
doubling, indicating that the pulsation of RR Lyrae 
stars may be even more complex. 

Our radial nonlinear hydrodynamic calculations 
demonstrated that the period-doubled fundamental 
mode may lose its stability with respect to the hrst radial 
overtone. Then it is straightforward to expect similar 
nonlinear dcstabilization with respect to low order nonra- 
dial modes as well. This process might provide the back- 
ground for the fact that additional pulsation frequencies 
occur more frequently in stars where period doubling is 
also detected. First overtone signals were identified in 
three other stars in the Kepler sample by IBenko et al.1 
(|2010f ). and the case of RRLyr suggests that more may be 
discovered. The period ratios of peaks that can be plau- 
sibly associated with the hrst overtone are for V354 Lyr 
(KIC 6183128): / x // = 0.729, and for V360 Lyr (KIC 
9697825): / x // = 0.721. For V445 Lyr (KIC 6186029), 
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IBenko et al.l ((2010) noted a very complex pattern of 
peaks while the detailed analysis of iGuggenberger et al.l 
(|2012f) concluded that among the additional frequen- 
cies, two can be found at ratios Jn/Io — 0.703 and 
/o//i = 0.730. Even more sta rs show the signatur e of 
the second overtone for which iMolnar et al.l (|2012h al- 
ready proposed the alternate idea of a three-mode reso- 
nance, where the equation 3/o + = /g may describe 
the relation between the three frequencies. 

The results indicate that interactions between radial 
and nonradial modes may i ndeed play a more crucia l 
role in the pulsation of stars. iBuchler fc Kollathl (|2011l ) 
showed that the resonance responsible for period dou- 
bling can cause Blazhko-like modulation in amplitude 
equations. Nonlinear asteroseismology promises to un- 
derstand these new findings, and might even help to un- 
lock the mystery of the Blazhko effect itself. 
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